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We examine the full counting statistics of electron transport through double quantum dots coupled 
in series, with particular attention being paid to the unique features originating from the level 
renormalization. It is clearly illustrated that the energy renormalization gives rise to a dynamical 
charge blockade mechanism, which eventually results in a super-Poissonian noise. Coupling of the 
double dots to an external heat bath leads to dephasing and relaxation mechanisms, which are 
demonstrated to suppress the noise in a unique way. 
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I. INTRODUCTION 

The manifestation of quantum coherence in finite sys- 
tems is the foundations of mesoscopic physics. Double 
quantum dots [l], Q , due to their inherent quantum co- 
herence, are widely accepted as promising candidates for 
building scalable qubits and quantum states detec- 
tors la, [7[ towards the realization of quantum computa- 
tion [8] . A great deal of effort has been invested to coher- 
ently measure, characterize and manipulate the quantum 
states in a double dot structure via coupling to a classic 
field or external surrounding 

It is well-known that a quantum system loses coher- 
ence due to coupling with a noisy environment. The in- 
volving issues have been the subject of intense research 
for many years [l6l - fl8j . Yet there is another important 
consequence of the system-environment coupling which 
renormalizes the internal energy of a quantum system. 
It has been revealed that, in a spin valve structure, the 
energy renormalization provides as an effective exchange 
magnetic field which leads to spin precession (l9l - [2l| . For 
a quantum dot Aharonov-Bohm interferometer, the level 
renormalization gives rise to additional dephasing of the 
quantum state |22| , as well as bias dependent phase shift 
and asymmetric interference patterns |23| . Even in solid- 
state quantum state measurement, the level renormaliza- 
tion was shown to play essential roles, and influence the 
measurement effectiveness crucially [2J, [25] ■ It is there- 
fore of vital importance to take this feature into account 
in order to correctly understand and analyze the electron 
transport properties. 

To have a specific example, we will investigate the en- 
ergy renormalization of electron transport through a se- 
rial double quantum dot system, as schematically shown 
in Fig.[TJ The analysis is based on the full counting statis- 
tics (FCS), which is capable of characterizing the cor- 
relations between charge transport events of all orders 
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FIG. 1: Schematic illustration for transport through serial 
double quantum dots. 



[26l |27I |. It thus serves as an essential tool superior to 
the average current in distinguishing various transport 
mechanisms involved [28|, |29(. We demonstrate unam- 
biguously that the level renormalization gives rise to a 
dynamical charge blockade mechanism, which eventually 
results in a pronounced super-Poissonian noise. Close 
relation between the super-Poissonian noise and the neg- 
ative differential conductance (NDC) owing to the level 
shift are revealed. It is further illustrated that the noise 
can be strongly suppressed due to coupling with an ex- 
ternal phonon bath. 

The paper is structured as follows. In Sec.|TT]we present 
the model Hamiltonian for the double dot transport sys- 
tem. Section lTm is devoted to the theory of FCS. The 
bias voltage dependence of the level renormalization and 
its influence on the FCS are discussed in Sec. lIVl which 
is then followed by the summary in Sec.lVl 



II. MODEL DESCRIPTION 

The system under study is schematically shown in 
Fig.Q] where the double dots are coupled in series. The 
entire system is described by the Hamiltonian 
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H = H B + H S + H'. 
The first part H B = £ Q=L R £ fc , CT £ 



(i) 

Cn,ka models 



the noninteracting electrons in the electrodes. Here c a ka 
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( c afccr) denotes the annihilation (creation) operator for 
electrons in the left (a= L) or right (a= R) electrode with 
spin a =t or I. The electron reservoirs are assumed to 
be in local equilibrium, so they are characterized by the 
Fermi functions /l/r(w). Electron transport is achieved 
by applying a bias voltage V, which is dropped symmet- 
rically at the left and right tunnel junctions. The bias 
then is modeled by different chemical potentials in the 
left and right electrodes Pl/r = ±V/2. Throughout this 
work, we set e = H = 1 for electron charge and the Planck 
constant, unless stated otherwise. 

The second part of the Hamiltonian depicts the cou- 
pled dots: 
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+ ^E( d LA- + 4,A.). 
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di„d n rr and n n 
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are the occu- 



pation number operators for dot a=L or R, with d atJ 
(d^ aa ) being the annihilation (creation) operator of an 
electron in the dot a with spin a. Each quantum dot 
is assumed to have only one spin-degenerate energy level 
-El/r within the bias window. One can parametrize the 
levels by their average energy E = (_El + En) / 2 and their 
difference e = E\, — Er, such that E h / n = E ± |e. VI 
accounts for the interdot tunneling between the coupled 
dots. Simultaneous occupation of one electron in each dot 
is associated with the interdot charge energy U'. Dou- 
ble occupation on the same dot cost the intradot charge 
energy Uq, which is assumed to be much larger than the 
bias voltage, such that charge states with three or more 
electrons in the double dots are prohibited. 

The third part H' = Y, a ,M<x<r d aa + d aafla) describes 
electron tunneling between the dots and electrodes, with 
f aa = ^ k t a kC a ka- The effects of the stochastic elec- 
tron reservoirs are encapsulated in the correlation func- 
tions C«(t -r) = (flMU'ir)) and C^t -r) = 
(f a a{t)f a<y ,(T)). Here, ((•■■)) = Tr[(- • • )p B ] represents 
the thermal average, with ps the local thermal equilib- 
rium reservoir state. Due to the serial geometry, electrons 
in the left (right) electrode can only transfer to the left 
(right) dot. The tunnel coupling strength of electrode 
a=L or R to the corresponding dot is characterized by 
the intrinsic line width r a (u;) = 27r^ fc \t a k\ 2 S(e a k — lj). 
In what follows, we consider only spin conserving tun- 
neling processes and assume flat bands in the electrodes, 
which yields energy independent couplings T a . 



III. FULL COUNTING STATISTICS 

The dynamics of the reduced system is described by 
the reduced density matrix p(t), which is obtained from 
the density matrix of the entire system by integrating 



out the reservoir degrees of freedom. Under the second- 
order Born-Markov approximation, it satisfies the quan- 
tum master equation [3Cj 



P (t) = -i£p(t) - np(t), 



(3) 



where the first part represents the internal dynamics on 
the double dots, with £(■••) = [Hs, (•••)]■ The second 
term accounts for the tunnel coupling between double 
dots and the external electrodes. Its detailed structure 
will be specified soon. The quantum master equation 
([3]) fully captures the dynamics of the reduced system; 
however, it is not adequate to describe the output char- 
acteristics. 

We unravel the reduced density matrix p(t) into com- 
ponents p( n ) (t) , in which "n" denotes the number of elec- 
trons that have been transferred to the right electrode. 
The resultant particle-number-resolved master equation 
reads [aMl, 



cr Iq=L,R 



^+p^A(+Jdl 



Rer 



H.c. 



(4) 



where A 



(±) 



(±) and A (±) 



iD aaa , (±£)}d aa i . Here, the involving spectral functions 
are defined as the Fourier transform of the bath correla- 
tion functions, i.e., 



C£U±£)= I dte ±ia c£l(t). (5) 



The dispersion function D^AzLC) then is determined 
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with V denoting the principal part. The spectral func- 
tions are associated with particle transfer processes, and 
the dispersion functions account for the coupling-induced 
energy renormalization of the dot levels. The latter has 
been neglected in previous work [36l - l40| . in which the 
Fermi energies of the electrodes are assumed to be far 
away from the electronic states of the dots. Later, it will 
be shown that the level renormalization can give rise to 
intriguing and important features in the output charac- 
teristics. 

By summing up all the components p( n ' (t) in Eq. (fj| , 
one recovers the unconditional quantum master equation. 
Consequently, the dissipative term in Eq. ((3]) is specified 
explicitly. The unique advantage of the particle-number- 
resolved master equation is its capability of establishing a 
close link between the reduced dynamics and the output 
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characteristics. By utilizing the conditional master equa- 
tion ((U), FCS characteristics can be readily determined, 
which enable us to get access to the complete information 
of transport. 

Let us start with the particle-number-resolved re- 
duced density matrix, which is directly related to the 
the probability distribution P(n, to) of having n electrons 
transferred through the system during the counting time 
to, i.e., P(n,to) = Tip^ n '(to), where the trace is over the 
reduced system states. The associated cumulant gener- 
ating function is defined as 

e -Hx) = Y J P{n,t Q )e-' mx , (7) 

n 

where \ is the so-called counting held. All cumulants of 
the current can be obtained from the generating function 
by performing derivatives with respect to the counting 
held 

{I k ) = -U-id X ) k Hx)\x=o- (8) 

The hrst four cumulants are related to the average cur- 
rent, the (zero-frequency) current shot noise, the skew- 
ness, and the kurtosis, respectively. 

To derive the cumulant generating function, we shall 
make use of the %-space counterpart of the number- 
resolved reduced density matrix g(\, t) = ^2 n p^ n \t)e mx . 
Its equation of motion, by employing the conditional 



master equation Q, reads formally 

g( X ) = C x g( X ), (9) 

where C x is totally determined by the dynamical struc- 
ture of Eq. Q . The formal solution can be readily 
obtained as g(x,to) = e Cxta g(x, 0). Straightforwardly, 
the cumulant generating function reads T{x,to) = 
— ln{Trp(x, to)}. Particularly, in the zero-frequency 
limit, i.e., the counting time to is much longer than the 
time of tunneling through the system, the cumulant gen- 
erating function is simplified to [27|, I4ll - l43 | 

^(x,*o) = -A roill (x)io, (10) 

where A m ; n (x) is the minimal eigenvalue of C x that sat- 
isfies A min | x ^ -> 0. 

IV. LEVEL RENORMALIZATION AND FCS 
ANALYSIS 

In the strong intradot Coulomb blockade regime, dou- 
ble occupation on the same dot is prohibited. The involv- 
ing states are restricted to: |0)-both dots empty, |L<r)- 
one electron in the left dot, |R<r)-one electron in the right 
dot, and |LcrR<7')-one electron in each dot, respectively. 
The quantum master equation ([3]) in this localized state 
representation reads 



Pha 



Po = - 
&(PrZ 

R<7 



2 ( r L + r R>o + r L (pLt + pu) + r R (pR t + pm), 
-p^-Fi 



PKa = m( P % - p&) - (r R + 2r+) PRa 

pLaRa' = — +r R )pLcrRcr' 



2r+)p L(T + r+po + r R (p Lff R<T + pl^r*), 

TrPO + T L (PLctRct + PhaRa), 



PRo' + f R pLo 



PRo = '^PRo + ^(PLa ~ PR.) - |(T£ + T R )p& - (T+ + T+> 
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(lie) 



r 



with spin a = {t?i} an( i & = — a - Here p s = (s\p\s) 
represents the diagonal element of the reduced density 
matrix. The off-diagonal elements p s s , = (s\p\s') de- 
scribes the so-called quantum "coherencies". The in- 
volving temperature-dependent tunneling rates are de- 
fined as T± = T a f^\E) and f± = T a f { a ±} (E + U% 
where /i^ (w) = {1 + e M^-^)/k B Tyi ig related to thc 
Fermi function of the electrode a— L or R. Here we are 
interested in thc regime A <C k^T (A = v 7 c 2 + 4f7 2 be- 
ing the eigenenergy separation), where the external cou- 
pling strongly modifies the internal dynamics, and the 
off-diagonal elements of the reduced density matrix have 
essential roles to play [H, |H, HiJ . The level separation 
is thus smeared by the temperature, and only excitation 
energies E and E + U' enters thc Fermi functions. 



By a close observation of the off-diagonal element of 
the reduced density matrix [see Eq. (|llcl) ]. it is found that 
the level detuning e is renormalized to 



e = e + A, 



(12) 



with A = Al — A R the energy renormalization. The level 
shift A a arising from tunnel coupling to the electrode a 
is given by 



K = <j> a (E) - 2ME + U') +4> a (E + U ), (13) 



with 
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2 27rfc B T J 



(14) 
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FIG. 2: (a) level renormalization A, (b) average current (I), (c) differential conductance G = d(I)/dV, (d) Fano factor 
F = (J 2 )/(J), (e) normalized skewness, and (f) normalized kurtosis versus the bias voltage for symmetric (Fl = Tr) and 
asymmetric (Fl = IOFr) tunnel couplings. Fhe bias is applied symmetrically at the left and right tunnel junctions, i.e., raise 
[il — V/2 and lower fiR = ~V/2 simultaneously. Each time when the chemical potential of the left electrode aligns with the 
energy needed for either single (E) or double occupation (E + U'), the level renormalization reaches its local extremum. Fhe 
plotting parameters are: e=0, r = T L + Tr . = 2fi, E = 10fc B T, U' = 20fc B T, and U = 100fc B T. 



Here \& is the digamma function. The involving intradot 
charging energy t/ serves as a natural cut-off for the 
energy renormalization [451 ] . 

The level renormalization is a genuine interaction ef- 
fect; it vanishes for Uq = V = 0. To clearly elucidate the 
effect of energy renormalization A, hereafter we assume 
e = 0. The numerical result of A is plotted in Fig.^a) 
as a function of bias voltage. Each time when the Fermi 
energy of the left electrode is resonant with the energy 
needed for single (E) or double occupation {E + {/'), the 
level shift reaches a local extremum. The energy renor- 
malization is sensitive to the tunnel-coupling asymmetry, 
i.e., A increases with rising Fl/F r ratio. 

The level renormalization gives rise to unique features 
in the transport current, as shown in Fig.^b). By ne- 
glecting the level shift, the current shows a typical step- 
like structure (see the solid and dashed curves). That 
is, each time when a new electronic level enters the bias 
window defined by the chemical potentials of the left and 
right electrodes, a current step occurs. The renormalized 
level detuning leads to suppression of the current (see 
the dotted and dash-dotted curves), particularly on the 
second current plateau, where at most one electron can 
reside on the double dots due to strong interdot and in- 
tradot charging energies. In this double-dot Coulomb 
blockade (DDCB) regime [H,[33|, an analytical expres- 
sion for the transport current can be obtained readily by 
utilizing Eq. ©, 



n 2 (2 + r R /2r L ) + r|/4 + A2' 



where the involving Fermi functions are approximated by 
either one or zero. Apparently, whenever the magnitude 
of the level detuning renormalization grows, the current 
reduces (the suppression of current close to E can be 
hardly resolved due to weak level renormalization and 
finite temperature) . The suppression of the current leads 
to regimes of NDC, as shown in Fig.^c). Actually, NDC 
has also been observed in many different contexts [4q . l47j . 

Double occupation on the system (one electron in each 
dot) becomes energetically allowed when the chemical po- 
tential of the left electrode crosses the excitation level 
E + U' . The current rises to the third plateau, which 
corresponds to the single-dot Coulomb blockade regime. 
The stationary current is given by 

(I) = 4(2r L + r R )r L r R » 2 

{ ' r L r R [(2r L + r R )2 + 4A 2 ] + 2(2r L + t k )*w ' 

(16) 

Again the current is suppressed whenever the level renor- 
malization grows [see also the dotted and dash-dotted 
curve in Fig.[2fb)]. 

Fet us now consider the second cumulant, which char- 
acterizes the width of the current distribution and is di- 
rectly related to the shot noise. Commonly, it can be 
expressed in terms of the so-called Fano factor F — 
(I 2 ) I (I) . The numerical result against the bias voltage is 
plotted in Fig.(2Jd) for symmetric and asymmetric tunnel 
couplings. At small bias (V <?C k-gT), the thermal noise 
dominates, and is described by the well-known hyper- 
bolic cotangent behavior which leads to a divergence of 
the Fano factor at V = 0. As bias increases but still well 
blow E, electron transport is exponentially suppressed. 
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FIG. 3: Schematic parallel two-level system obtained by a 
unitary transform of the double quantum dots coupled in se- 
rial as shown in Fig.[T] The tunneling amplitudes are effec- 
tively modulated by the level renormalization. 



Tunneling events are uncorrelated and the noise exhibits 
Poissonian statistics [28| . The transport through the sys- 
tem becomes energetically allowed when the bias is fur- 
ther increased to the DDCB regime. Owing to the level 
renormalization, the Fano factor exhibits clear enhance- 
ment at bias close to the excitation energies E and E+V. 
Hence the Fano factor proves to be much sensitive to the 
internal energy than the average current. 

Noticeably, a pronounced super-Poissonian noise is ob- 
served in the DDCB regime [see the dash-dotted curve in 
Fig.^d)]. By evaluating the minimal eigenvalue of £ x , 
one obtains from Eq. ([8]) the analytical expression of the 
Fano factor 



f = i-4iv n 



: r R (r R + 6r L r R + sn 2 ) + 4A 2 (r R - 2r L ) 
[2ft2(4r L + r R ) + r L (r R + 4A2)]2 ■ 

(17) 

Unambiguously, super-Poissonian noise is expected when 
the second term is negative. This is satisfied under the 
conditions of A ^ and T R < 2r L . In this case, the 
tunnel-coupling to the right electrode is not very strong, 
and thus the Coulomb interactions are more effective. 

In the DDCB regime, states with two or more elec- 
trons in the double dots are prohibited. We perform a 
unitary transformation to diagonalize the reduced sys- 
tem Hamiltonian as (spin indices are suppressed here in 
order to simplify the discussion) 



h s = |A(H-><+! - |-x-|), 



(18) 



where A = VA 2 + 4S1 2 , with the level renormalization 
being properly accounted for. The serial double dots 
are thus mapped onto a parallel two-level system, as 
schematically shown in Fig.[3J Here, the energy eigen- 
states are denned as 



|+) = sin -|R) + cos -|L), (bonding state) (19a) 
9 9 

|— ) = cos -|R) — sin-|L), (anti-bonding state) (19b) 



where 9 is introduced via sin 9 — 2^1/ A and cos 9 = A/ A. 
As a result, the system-electrode coupling Hamiltonian 



is recast to 



cos o c L + tnk sin -c\ k |0) (+| 



+ (t Rk cos ^c Rfe - t Lk sin |c£k) |0)(-|| +H.c. 

(20) 

Apparently, the effective amplitudes (as shown in Fig. [3]) 
of electron tunneling through the bonding and anti- 
bonding states are affected notably by the level renor- 
malization. 

For A = (9 = it/2) and symmetric tunneling ampli- 
tudes (thk = i R fc), the noise of electron transport through 
the bonding and anti-bonding states is suppressed due to 
the Pauli exclusion principle. It leads thus to a sub- 
Poissonian statistics, as shown by the solid curve in 
Fig. Eld). Finite tunnel-coupling asymmetry enhances 
the degree of correlation in transport. Yet, the noise 
can not exceed the Poisson value by increase the tunnel- 
coupling asymmetry alone [cf. Eq. ([T7| ]. Correlation of 
transport can be crucially increased by the level renor- 
malization. As the level shift A grows, electron transport 
through the bonding and anti-bonding states strongly 
modulate each other, and a dynamical channel blockade 
mechanism is developed [H], H-[5l| . It gives rise to the 
bunching of tunneling events, and eventually results in 
the super-Poissonian noise, as shown by the dash-dotted 
curve in Fig.[2](d). 

The occurrence of the dynamic charge blockade leads 
generally to the suppression of the current, which finally 
gives rise to the NDC. However, the NDC does not nec- 
essary imply the super-Poissonian noise, as shown by the 
dotted curves in Fig.JHc) and (d). This seems to be at 
variance with that in Ref. [52| . where electron trans- 
port through a multi-level quantum dot is investigated. 
There the NDC is found to be always accompanied with 
super-Poissonian noise. NDC was also observed in Ref. 
[17j . where the occurrence of NDC was associated with 
reduction of charge accumulation due to bias dependent 
tunneling rates. The NDC in the present case is a pure 
energy renormalization effect, as will be explained below. 

For the present double dot system, the rates of tun- 
neling through the bonding and anti-bonding states (as 
shown in Fig. [3J depend not only on the tunneling am- 
plitude (t a k), but also on the level renormalization A. In 
the case £l/c = iiifc, enhancement of the level renormaliza- 
tion results in a decrease of sin | . The rates of electron 
tunneling from left electrode to the anti-bonding state 
and that from bonding state to the right electrode are 
both suppressed. As a result, the current is reduced and 
NDC occurs [see the dotted curve in Fig. [2c)]. How- 
ever, the noise does not exceed the Poissonian value due 
to suppressed transport through both channels. Now, 
consider an increase of tuc, which enhances the rate of 
electron tunneling from the left electrode to the double 
dots. In the presence of a strong energy renormalization, 
the rate of electron tunneling from the bonding state to 
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FIG. 4: (A) Fano factor and (b) normalized skewness versus 
bias voltage for different dissipative couplings 77. The tun- 
nel couplings to the left and right electrodes are asymmetric 
(Fl = lOrrt) and the energy cutoff for the phonon bath spec- 
tral density is u c — 5T. The other parameters are the same 
as those in Fig. 



the right electrode remains very low. In the limit where 
the Coulomb interactions prevent a double occupancy 
of system, there is competition between the two trans- 
port channels. Consequently, the slow flowing of elec- 
trons through the bonding state modulates that through 
the anti-bonding state, which gives rise to a bunching 
of tunneling events and eventually leads to the super- 
Poissonian noise and NDC, as displayed by the dash- 
dotted curves in Fig.^c) and (d). 

We are now in a position to examine the third and 
fourth cumulants, which characterize respectively the 
asymmetry and sharpness of the current distribution. 
The normalized skewness is displayed in Fig.[2Je). In 
the DDCB regime, it exhibits notable enhancement due 
to energy renormalization. In particular, a pronounced 
super-Poissonian behavior is observed for strongly asym- 
metric tunnel couplings, as displayed by the dash-dotted 
curve in Fig.^e). The kurtosis, as shown in Fig.[2]T), can 
be either enhanced or reduced by the energy renormaliza- 
tion, depending sensitively on the tunnel coupling asym- 
metry. The peak width of kurtosis, in comparison with 
that of Fano factor and skewness, is reduced obviously; 
it thus reflects more precisely where the resonance and 
the extremum of energy renormalization are achieved. 

Finally, let us turn to the influence of external phonon 
bath which leads to relaxation and dephasing in the sys- 
tem 0, |53| . The corresponding rates are given respec- 



tively by |4J, |48|, |4£ 



and 



4- coth ( — — ! ± 1 
A \2k B T 



777T- 



(21a) 



7 P h = 2 7 r|jJ ph (A)coth(-^) 



?77T- 



-, (21b) 



where J p h(w) = rjuie^^'^ is the Ohmic spectral density 
of the heat bath, with the dimensionless parameter 77 
reflecting the strength of dissipation and ui c the Ohmic 
high energy cutoff. 

The calculated Fano factor and normalized skewness 
versus bias voltage are plotted in Fig. 2] for different dis- 
sipative couplings r\. Strong suppression of the cumu- 
lants is observed in the DDCB regime, where double oc- 
cupation on the system is prohibited. In the absence of 
electron-phonon coupling, electrons are transferred co- 
herently. The phonon bath coupling, on one hand, leads 
to dephasing of the quantum state, and electrons tend to 
be tunneled sequentially, which then causes the reduction 
of noise in the DDCB regime [48]. On the other hand, the 
involving phonon emission and absorption processes give 
rises to a relaxation mechanism, which further suppresses 
the noise. The noise suppression is particularly notable in 
the bias close to E+U', where the energy renormalization 
is prominent, as shown by the dotted curves in Fig.^ta) 
and (b) for a strong dissipative coupling 77 = 0.02. These 
features demonstrate that both Fano factor and skewness 
are sensitive tools to the phonon bath induced dephasing 
and relaxation. 



V. CONCLUSION 

In summary, we have investigated the full counting 
statistics of electron transport through double quantum 
dots coupled in series, and paid particular attention to 
the unique features arising from the level renormaliza- 
tion. It is found that a dynamical channel blockade mech- 
anism is developed purely due to the energy renormal- 
ization, which eventually leads to a pronounced super- 
Poissonian noise. Our results demonstrate unambigu- 
ously the importance of the level renormalization when 
investigating the transport properties of a double dot 
structure. Negative differential conductance due to level 
detuning renormalization is observed, and its relation 
with the super-Poissonian noise is revealed. Coupling 
of the double dots to an external phonon bath leads to 
dephasing and relaxation mechanisms, which are shown 
to suppress noise notably Furthermore, double dot sys- 
tems, as recently shown in Refs. [HI, (55[, are good can- 
didates for entanglement generation. In this context, the 
present full counting statistics have the potential to fa- 
cilitate the identification and characterization of entan- 
glement originating from different sources. 
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